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Abstract 

One-dimensional nonlinear lattices with harmonic long range interaction 
potentials (LRIP) of inverse power type are studied. For the nearest neigh- 
bour nonlinear interaction we shall consider the anharmonic potential of the 
Fermi-Pasta-Ulam problem and the (f>'^+(f>'^ as well. The continuum limit is ob- 
tained following the method used by Ishimori [1], and several Boussinesq and 
KdV-type equations are found with supplementary Hilbert transform terms. 
This nonlocal terms are introduced by the LRIP. For the (f>^ + (f>'^ nearest 
neighbour potential, the continuum approximation turns out to admit exact 
bilinearization in Hirota formalism. Exact rational nonsingular solutions are 
found. The connection with perturbed KdV equations is also discussed, and 
the dynamics of the KdV soliton is studied for one of these equations. 

Introduction 



Several years ago Ishimori [1] has discussed in the continuum approximation a one- 
dimensional nonlinear lattice with a long range interaction potential of Lennard- 
Jones type. This was done expanding the lattice displacements in Fourier series, 
performing the summations in the real space, and finally taking into account only 
the first terms in the wave vector expansion - the long wave limit. His conclusion 
was that the long range interaction manifests mostly in the linear dispersive terms, 
the effect on the nonlinear terms being less dramatic. Several integro-differential 
equations were obtained, some of them containing an Hilbert transform. 

From another point of view, in the context of wave propagation in shallow water, 
Whitham [2] has extended the well known Korteweg-de Vries equation by including 
a nonlocal (integral) dispersive term. An example from plasma physics is given by 
the nonlinear ion acoustic waves where a Hilbert transform term is introduced into 
the KdV equation to describe the Landau damping [3]. 

Recently the interest in the study of nonlinear evolution equations with non- 
local dispersive terms was stimulated by the discovery that their solutions have 
very interesting picking and breaking phenomena [4]. On the other hand it is well 
known that several integro-differential equations of Benjamin-Ono (BO) type [5] 
are completely integrablc and the main feature of their integrability character is the 
existence of rational solutions (algebraic solitons). 

Following Ishimori [1] we shall consider a one-dimensional lattice with long range 
harmonic interaction between the atoms, the nonlinear interaction being restricted 
between nearest neighbours. The long range harmonic interaction potential is writ- 
ten 
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<2^m,n'^m'n{Um Un) , (1) 

where ^ ^ ^ 

2 C,{p) \m — n|P 

and in the summation the term m = n is excluded. Here u„, are the dimension- 
less lattice displacements at the lattice points m,n, p >2 is an integer describing 
the spatial decreasing of the long range interaction, and J measures its strength. 
In the definition of Jmn we introduced the Riemann sum 



oo ^ 
m=l 

and consequently 

T:^Jmn = J. (4) 

Also a long range potential with alternating sign will be considered. For this 
Jmn is given by 

2 |m — n|P 

and in this case 

Y^^Jmn = J(l - (4a) 
The nonlinear interaction is assumed only between nearest neighbours 

^ - u„) (5) 

and the simplest form for V is 

V{u) = ^^u\ (6a) 

a being a parameter describing the strength of the anharmonic term. The expression 
(6) is the same as in the famous Fermi-Ulam-Pasta problem. We shall discuss, also 
the following form 

Viu) = \u- + + ^u' (66) 

which is an anharmonic potential introduced by M. Wadati [11]. Obviously the 

continuum approximation for the equations of motion will involve Boussinesq type 
equations with long range interaction corrections. 

Due to the inverse power form of the long range kernel, continuum limits of the 
long range interaction terms turn out to be in the form of a Hilbert transform and 
its derivatives [1]. Also the continuum limits of the local terms, restricted at the 
waves in only one direction, leads to the appearance of the KdV [7] equation for the 
local potential (6a), and KdV+mKdV [11] equation for the local potential (6b). In 
the first case, assuming J to be small a perturbatively approach is done in order to 
find what happened with the 1-soliton solution of KdV type. It is found that for 
p = 2 the first 5 conservation laws for KdV 1-soliton remain unchanged in the first 
order of perturbation. 

In the second case the presence of the Hilbert transform which is a nonlocal 
dispersive term leads us to the possibility to consider rational solutions which ap- 
pear in nonlinear integro-differential equations such as BO-hierarchy [5] and nonlo- 
cal AKNS-hierarchy [5]. These rational solutions behave like solitons and interact 
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elastically with no phase shift. For local (1+1) dimensional nonlinear evolution 
equations both continuous and discrete the rational solutions also exist but their 
character is rather singular [12], [13], [14], [17]. They have no free parameters and 
there are very few cases in which thciv are not singular [12], [15], [16] and, moreover, 
they interact elastically with no phase shift [16]. 

It turns out that for p = 4 and (6b) the continuum limit of the equation of 
motion admits exact bilinearization in Hirota formalism for certain values of the 
parameters a and (3. The corresponding rational solutions are exactly the same as 
for mKdV with nonvanishing boundary conditions [12], [15]. In this way we are 
expecting that this system to be completely integrablc although it does not belong 
to the BO-hierarchy or Nonlocal AKNS-hierarchy either. 

Continuum approximation for the equations of motion 



The equation of motion writes 

= V (U„+1 -Un)-V {Un - Un-l) - 2'E^Jmn{Un - Wm)- (7) 

In the linear approximation it becomes 
Looking for solutions of the form 



the following dispersion relation results 

a;2(fc)=4sin2^ + 2J^Fp(fc), (9) 

where 

„ , 1 — cos mk 

= E — — • (10) 

m— 1 

This is an even function of k and can be calculated for each value of p [6] . In the 
following the results for several values of p will be presented. 

Ccise p=2 



F2{k) = \\k\-\k\ C(2) = y. (11) 



Then 



a;^(fc)=4sin^ ^ + 2J - ^fc^ , ^ (12) 



which in the long wave limit is approximated by 



^■^{k) ~ + (1 _ _ xk^ + 0{k'), (13) 

TT TT 
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where A = ^ is a measure of the dispersion. This dispersion relation corresponds 
to a hnear partial differential equation, namely 

Utt = aH{ux) + c^Uxx + ^u^xxx, (14) 

where the subscript indicates the derivative with respect to the corresponding vari- 
able, and the following notations were introduced 

a= — , = 1 5-. (15) 

Here H{f{x)) is the Hilbert transform of the hmction f(x) 

H{f{x)) = -P r°^ 4^dx, (16) 

TT X - X 

and we have used the relation 

Hie''''') = isgnk e''^^. (16a) 
The contribution of the anharmonic term av? in the equation of motion is 

-a(u„+l - Unf + a{Un - Un-lf = -a(u„+i - U„_i)(m„+i - 2Un + Un-l) 

and in the continuum limit this goes to 

-2aUxU^x = -a{ul)x. (17) 
Then the complete nonlinear partial differential equation in the continuum limit is 
Utt - C^Uxx + Ot{ul)x - ^Uxxxx - aH{ux) = 0. (18) 

This is a Boussinesq type equation with a correction term in the form of a Hilbert 
transform because of the long range harmonic interaction. 

Case p=3 

In this case (see the Appendix) 

Fsik) = In |fc| + ^-k' ^k^ + 0{k% (19) 
The dispersion relation becomes 

and the corresponding Boussinesq type equation is 

Utt - c^Uxx + a{ul)x - >^Uxxxx + aT{uxxx) = 0, (21) 

where 



c2 = l + 



3J 
2C(3) 



J 



and T{f{x)) is the following integral operator [1] 

T{f{x)) = ^ p sgn {x - x) (in \x - x| + 7) f{x)dx (23) 
7 being the Euler constant. In deriving the above equation we have used 



r(e''=^) = ililMe*'^^ (23a) 



Case p=4 

2-1 4 

F4(fc) = — fc2_^|fc3|^^^4 C(4) = — , (24) 

^ ^ 12 12' ' 48 ' ^ 90' ^ ' 

and the following equation is finally obtained. 

Utt - C^Uxx + a(Wx)x - AUxxxx - aH{Ua;xx) = 0. (25) 

Here 

2 , 15J 15J , 1 45J, , 



For an alternating in sign long range interaction potential, when J^n is given 
by (2a), the same technique gives slightly different results. Instead of (10) Fp{k) is 
now defined as 



f;^"Hfe) = E(-ir"^ (26) 

m — l 

As an example we shall consider the case p = 2 when [6] 

Fj'^)(fc) = ^P. (27) 

Consequently the dispersion relation in the long wave limit takes a very simple form 

ij^{k) = c^k'^ - Afc* + 0(fc*^) (28) 
= 1 + 3 J 

and the corresponding nonlinear evolution equation in this continuum limit becomes 

Utt - c^Uxx + ^Uxxxx + a{ul)r, = 0. (29) 

No non-local corrections appear and the long range interaction has a small influence 
modifying only the wave velocity. 

It is well known since the paper of Zabusky and Kniskal [7] that the KdV 
eq. is obtained as the uni-directional long time approximation of the Boussinesq 
equation. We have to assume that there is a competition between the dispersion, 
nonlincarity and wc have to add now the new term in the form of the Hilbcrt 
transform, competition which manifests at long times. In order to describe it we 
shall assume 

a — > ea, A —>■ eX, a — > ea, (30) 
where e is a small parameter. The solution u{x, t) is written as 
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u{x,t) = f{^,T) + eu^^\x,t) (31) 

where 

^ = X — ct 

T = et (32) 
are the "slow variables". Then in the first order in e from (18) we get 

4'^ - = 2c/,. - a(/2), + A/,,,, + aH{f^ ) . (33) 

As the right hand side depends only on ^ and not on i] = x + ct it has to vanish, in 
order to prevent the linear rising in r] of u^^\x,t). Defining 



9= 



r=l. (34) 

we finally obtain 

qr - 6qq^ + Xq^^^ + aH{q) = o (35) 

which in a KdV equation with an additional term (long range correction) in the 
form of a Hilbert transform. 

If the same procedure is applied to the equation (25) we get 

Qr - 6qq^ + Xq^^^ - aH{q^^ ) = o. (36) 

Here the correction term has the same form as the Hilbert transform appearing in 
the Benjamin-Ono equation [5], which is actually obtained if the dispersive term is 
dropped out. 

A quite surprising result emerges in the case of long range interaction with 
alternating signs. As shown before in this case the continuum limit is the usual 
Boussinesq equation (29) and from it the KdV eq. is immediately derived. It seems 
that the alternating signs in the long range interaction potential (we can call it an 
" antifcrromagnetic" interaction) has a very small effect on the nature of elementary 
excitations in the long wave limit. This conjecture has to be verified for other values 
of p in (2a). 



Perturbatively approach 



Now for small values of the parameter "a" one can ask ourselves in what way the 
soliton of the KdV eq. is perturbed. We shall discuss briefly this problem for 
eq. (35) (we shall put A = 1) following one of the methods used to the analytical 
description of soliton dynamics in nearly integrable systems [8]. The approach is 
based on calculating the changes ( evolution in time ) of the first conservation laws. 
It is easily seen from eq. (35) that the first two, the mass and impulse, remain 
unchanged 

/+(X) 
q{^,T)d^ = ct. 
-OO 

/ + 00 
q^{^,T)d^ = ct. (37) 
-OO 
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In order to find the time evolution of the KdV soliton energy we have to multiply 
eq. (35) by 



„ 2 9 
3a 



and integrate over ^. One obtains 



^ + a da^' + 4?)H{q{0) = (38) 



dr J_oo * 9^ 

where for the one soliton solution 

gfc(^r) = -2k^sech^ {k{^ - 4A;V)) 



Ek= 1^°^ (U^+Q')d^ = -'^k^ (39) 



.2 « 

The second term in (38) separates in two integrals, namely: 



7,=ipr r,«*'i^ ,40) 



and 



After a partial integration on ^ in /2 we get 



h = —P / di (41a) 



But for the KdV soliton 
Then using the fact that 



+ 3g^ (42) 



pp^^'i(|)j|)=„. (43) 

J — oo J — oo ^ > 

relation which was already used for proving the momentum conservation (37), we 
obtain 

/i + /2 = (44) 

and consequently also the soliton energy is conserved in a first order perturbation 
calculation. 

We can perform the same procedure for the following conserved quantity: 

Ti = 5/ + lOggf + q^^ 

Multiply eq (35) by: 

and integrating over ^ using the same procedure it is found that T4 remains con- 
served. Lenghty but straightforward calculations lead to the same result for the 
next conserved quantity: 

n = 219^ + 1059^91 + 21991^ + l^lii 
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We can proceed to any conserved quantity but the calculations are very much 
involved. 

This is quite surprising because for an equation not very different from (35), 
describing the Landau damping of an ion acoustic wave [3], namely 

Uj- — 6nn^ + n^^j + aH{n^) = (45) 

only the first conservation law (particle conservation) is not destroyed, the other 
ones decaying in time. 

In the same time wc have to stress that our result is valid only in the first 
order of a perturbation theory and the resulting conclusions have to be understood 
only in this sense. More elaborated techniques have to be used to go beyond this 
approximation [8], [9]. 

That the situation is much more complicated was proven ten years ago by Birnir 
[10]. Discussing an equation of the same form as (35), with the coefficient a in front 
of the Hilbert transform a periodic function, satisfying mild conditions, he was able 
to prove the existence of chaotic rational solutions. They look like their integrable 
brothers, but are moving in a chaotic fashion and also can disappear and reappear 
[10]. This result shows how complex the problem of the perturbed KdV equation 
(35) is. 



Exact bilinearization 



In this case the local potential is given by: 

V{u) = + ^u' + {u^ (46) 

where wc assiimc /? > For p = 4 the continuum approximation will have one more 
nonlinear term due to the anharmonic term in the potential, so: 

wtt - ^u^x - a{ul)x - P{ul)x - Amxxxx - aHu^xx = (47) 
Following the same procedure, we consider: 

u{x,t) = f{^,T) + eu^'\x,t) (48) 

and one gets: 

u[l^ - c^uW = 2c/jT + a(/|)j + /3(/|)i + Xf^a + aif/«« = (48) 
For T = T/2c the equation becomes, after one integration: 

fr + afl+pfl + Xf^^^+aHf^^ = (49a) 

which is a potential version of KdV+mKdV equation with nonlocal dispersion. For 

the sake of simplicity wc shall consider J such that X = a = 1. One can see for 
r = T/2c and g(C,r) = i/j, (48) becomes: 

qr + 6aqq^ + 27f3q'^q^ + + Hq^^^ = (496) 
In order to bilinearize the equation (49a) we shall use the following substitution: 
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where 



N 



g.{^,r) = ]l{^-Zn{T)) (51a) 

n=l 

N 

<?+(C,r) = n(e-<(r)) (516) 



71=1 

and Imznir) > 
With this ansatz 

^fii = log (5+5-) 
and introducing in (49a) the following expression appears: 

glg'^it'Dr - D| + iD|)5+ . g_] -{a- l)glg^_{T>^g+ • g.f- 

-3iglg^_(D^g+ • 5-)(D|5+ o <?-) + *(2 - /?)5+5-(D«5+ • 9-f = (52) 

where D"/ • g = (9a; — dy)"f{x)g{y) \x=y are the Hirota bilinear operators. 
If a = 1, /3 = 2 and decouples (52) one gets: 

{iBr - D| + iD|)5+ •g-=0 
D|5+.5_ =0 

which relaxes to 

(D^ + D|)5+.5_ =0 
D|ff+.5_ =0 

and it is nothing but bilinear form for mKdV equation. There is no solution of the 

form (51a), (51b). 

If /3 = 2 and 1 one obtains: 

{iBr - D| + ir>l)g+ •g-=0 



B'^ - -{a - ] g+ • g- = (53) 



The equation (53) is equivalent with 



iT>r - ^(a - 1)D4 +i-Dl]g+*g_=0 



D2--(a-l)Djj5+.5_ = (53) 

which represents the bilinear form for the following mKdV equation with non- 
vanishing boundary condition: 

Ur + 3(q! — l)u^U^ + M^jj = (54) 

and u — > —1/3 when ^ ^ ±(X). This equation is completely integrable and it admits 
A/'-soliton solution. By the long-wave limit procedure of Ablowitz and Satsuma it 
admits also nonsingular real rational solutions in the form (51a) and (51b) and 
accordingly they are solutions of (49). Due to the fact that (54) admits rational 
nonsingular soluitons for every odd N, we can conclude that the set of rational 
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solutions for (49) is the same with the set of rational solutions for (54). The 1- 
rational and 3-rational solutions for (49a) are given by: 



3i 



a-1 



+ d 



and 



27 ^ 9^ 



(a - 1)2 



a-1 
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{a-iy 



+ d 



where 0{^,t) = + {a — + c) and c, d are arbitrary constants. Thus, 

q = i(9^ log(g-|_/(7_) is a nonsingular rational solution We can see that only for 
A'' = 1 the solution is a solitary wave. The velocity is fixed and the 3-rational 
solution is not a superposition of solitary waves. 



Appendix 



In calculating Fk we start from 



and the integral 



sin(mfc) Z"*^ , . t . , 

^ - ' ln(2sin-)dt. 



sin(mfc) f 
^1 ~ Jo 

m=l " 



We have 



But 



and 



F^(k) = - f dk f \n{2smUdt 
Jo Jo 2 

f'' t 
= - / (A;-i)ln(2sin-)dt 
Jo 2 

= _^ln(2sin^) + l^ [kt -*^) cot ^dt. 

ln(2sin-) ~ Infc + ln(l - — + ...) ~ Infc - — + O(fc^) 
2 24 24 



Then, by straightforward integration, we get {k > 0) 

Fs{k) = -^lnk+le + ^k' + 0{k% 
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